Level and Eigenfunction Statistics in Billiards with Surface Scattering 
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Statistical properties of billiards with diffusive boundary scattering are investigated by means of 
the supersymmetric cr-model in a formulation appropriate for chaotic ballistic systems. We study 
level statistics, parametric level statistics, and properties of electron wavefunctions. In the universal 
regime, our results reproduce conclusions of the random matrix theory, while beyond this regime we 
obtain a variety of system-specific results determined by the classical dynamics in the billiard. Most 
notably, we find that level correlations do not vanish at arbitrary separation between energy levels, 
or if measured at arbitrarily large difference of magnetic fields. Saturation of the level number 
variance indicates strong rigidity of the spectrum. To study spatial correlations of wavefunction 
amplitudes, we reanalyze and refine derivation of the ballistic version of the cr-model. This allows 
us to obtain a proper matching of universal short-scale correlations with system-specific ones. 



PACS numbers: 73.23.Ps, 05.45. +b, 73.20.Dx 



I. INTRODUCTION 

Chaotic cavities, commonly understood as quantum 
systems whose classical analogues exhibit chaotic dy- 
namics, have become a common object of research in 
condensed matter physics. Experimentally, features of 
chaotic motion appear, among others, in quanturn-,dotsEJ, 
microwave cavitiesa, and small metallic clusters^. For 
defmitcness, we talk below about electrons in quantum 
dots (billiards). 

On the theoretical side, properties of chaotic cavities 
are subdivided into universal and non-universal. By uni- 
versal we mean physical quantities which only depend 
on the global symmetry of the system (like the time re- 
versibility) and possibly trivially scale with the size of 
the system, but are not influenced by any cavity-specific 
details of dynamics of the electron motion. At the same 
time, these quantities differ for systems with chaotic and 
integrable classical analogues, and are therefore concep- 
tually important when one discusses signatures of chaotic 
behavior in quantum systems. Examples of universal ef- 
fects include low-frequency level statistics, leading order 
conductance and shot noise. These properties are by now 
well understood and described by various types of ran- 
dom matrix theory (Gaussian ensembles of Hamiltonians 
for closed or almost closed systems or circular ensembles 
of scattering matrices for open systems). 

On the other hand, non-universal quantities, such 
as, for instance, level correlation function for high en- 
ergy separation or eigenfunction correlations at large dis- 
tances, are determined by sample-specific details of elec- 
tron motion. These non-universal quantities thus dis- 
criminate between the behavior of individual billiards. 

A standard tool for treating fluctuations of the density 
of states (DOS) in chaotic systems is the real space path 



integral approach. Within this method, the DOS corj: 
lation function is given by Gutzwiller's trace formula 
which has a form of a sum over periodic orbits of a specific 
billiard. To evaluate the formula explicitly in the non- 
universal regime, one has to resort to numerical treat- 
ment. 

In this paper we will use an alternative approach 
which has attracted considerable interest recently, the 
ballistic cr-model. It generalizes the supersymmetric cr- 
model, whici|-|proved to be very successful lor disor- 
dered metalsQl], to ballistic disordered systemsBfl. In the 
framework of this approach, all non-universal quantities 
are expressed through eigenvalues and eigenfunctions of 
the Liouville operator, which introduces non-universality. 
It has also been conjectured that the same cr-model in the 
limit of vanishing disorder describes statistical properties 
of spectra of individual classically chaotic system. This 
conjecture was further developed in Refs. where 
the cr-model was obtained by means of energy averaging, 
and the Liouville operator was replaced by its regular- 
ization — the Perron- Frobenius operator. The progress 
along this direction is complicated by the fact that the 
eigenvalues of the Perron-Frobenius operator for many 
systems are unknown, while its eigenfunctions can be ex- 
tremely singular. 

We thus conclude that it is highly desirable to have an 
example of a chaotic system which is analytically solv- 
able. We would expect that universal properties of such 
an example will conform with the predictions of the ran- 
dom matrix theory (RMT), whereas explicit expressions 
for non-universal quantities would improve our under- 
standing of properties of chaotic cavities. 

Currently, we are unaware of such an example with 
chaotic dynamics. However, one can instead treat sys- 
tems with surface disorder, which leads to diffusive scat- 
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tering at the boundary of a billiard. This model mimics 
the behavior of a system in the hard chaos regime: as a 
result of surface disorder any two arbitrary close trajecto- 
ries spread apart after the first collision with the surface. 
This must be contrasted with slightly distorted integrable 
billiard with the typical spatial scale of this distortion 
being of the same order as the size of the systemllialla. 
Those systems, termed by the authors rough billiards, 
exhibit slow diffusion over angular momentum. Systems 
with surface disorder are also different from intertable 
systems with bulk disorder in the ballistic regimeLTlI^ . 

Studies of level and eigenfunction statistics00 have 
shown that, indeed, universal properties of a billiard with 
diffusive surface scattering agree with the RMT predic- 
tions. At the same time, non-universal features have been 
found which reflect the classical ballistic dynamics in the 
billiard, j— Subsequently, analytical results for persistent 
currentsE3 and transport properties^ of chaotic cavities 
have been obtained in the same model. Numerically, bil- 

in connec- 



statistics, 

Refs. Il-lll 
cites having 
which mod- 



liards with surfacerdisarder have been stui 
tion with the levele-ffia and eigenfunctionE 
and magnetoconductancecH. To this end 
consider a lattice model with the boundar 
random energies. A treatment of Ref. |28| 
els surface disorder by cutting off the boundary cites by 
confining potential, and eventually proceeds with numer- 
ical evaluation of persistent currents, seems to describe a 
similar physical situation. 

In this article, we perform a systematic analytical 
study of a circular billiard with boundary disorder, based 
on the cr-model approach. The paper is organized as 
follows. Section |l| presents the cr-model for a circular 
billiard with diffusive boundary scattering. We subse- 
quently use this approach to derive the results for level 
statistics (Section QI), parametric level statistics describ- 
ing variation of individual energy levels with applied 
magnetic field (Section [Tv| ), and eigenfunction correla- 
tions (Section [y]). In Section [Vl| we generalize the prob- 
lem, imposing the mixed boundary condition, instead of 
purely diffusive reflection. This boundary condition en- 
ables us to model a broader class of chaotic systems, 
where the lowest Lyapunov exponent is parametrically 
different form the inverse time of fl ight . The obtained 
results are summarized in Section |VTj , where we also 
present a discussion of some op en p roblems. A brief ac- 
count of the results of Sections |l| and |v| has been pre- 
viously given in Ref. ^lj; the level statistics (Section III) 
were independently studied in Ref. 



II. CIRCULAR BILLIARD WITH DIFFUSIVE 
BOUNDARY SCATTERING: a-MODEL 
APPROACH 



A. General considerations 

We consider a 2D circular billiard of a radius R, which 
is clean (ballistic) inside, and contains some disorder (to 
be specified below) at the boundary. Our starting; point 
is the cr-model for ballistic disordered systemsOfi. The 
effective action for this model has the form 



F[g(r,n)} = - 



drSti 



iu)A(g(r)) 



1 



2r(r) 



2v F (AT- 1 nVT)] 



(1) 



Here the semi-classical Green's function g(r,n) inte- 
grated over energies is a 4 x 4 supermatrix which de- 
pends on the coordinate r and direction of the momen- 
tum n. To simplify the presentation we consider the case 
when the time-reversal symmetry is broken in the quan- 
tum problem but is preserved in the classical one, which 
can be achieved e.g. by ap plyin g a weak magnetic field. 
(See discussion in Section |IVC ). The angular brackets 
denote averaging over n: (0(n)) = J dnO(n) with the 
normalization J dn = 1, and the supertrace is defined 
as trace of boson-boson block minus trace of fermion- 
fermion one. The matrix g is constrained by the condi- 
tion g(r,n) 2 = 1, and generally can be represented as 
g = TAT -1 , with the matrix A = diag(l, 1, -1, -1) dis- 
criminating between retarded and advanced components 
of the Green's function. As usual, vf and v — m/(2Tr) 
denote the Fermi velocity and the density of states at the 
Fermi surface, respectively; r(r) is the (position depen- 
dent) elastic scattering time, which originates from the 
disorder. We use the units with ft = 1 in the rest of the 
paper. 

When all the disorder is at the boundary the scatter- 
ing time t(t) must be chosen in a way that it is infinite 
everywhere except for a thin layer around the boundary. 
In the consideration below this term only modifies the 
boundary condition. 

The action £j]) differs from that in the cr-model for dif- 
fusive systemsETin two respects: First, the Green's func- 
tion g is defined in phase space; second, Eq. ([!]) is lin- 
ear in gradients, whereas its diffusive counterpart is of 
second order in spatial derivatives. Despite this differ- 
ence, methods developed for the calculation of level and 
eigenfunction statistics in diffusive systems can be ap- 
plied here. Indeed, these properties are governed by the 
structure of the action in tha.vicinity of the homogeneous 
configuration of the g-fieldEj, g(r,n) — A (see Refs. 
for review). In this case the action may be considerably 
simplified. Writing T = 1 — W/2 + . . ., we find the action 
to leading order in W, 



F [W] 



~~2 



drdnStr 



W21 (k - iu\ W 12 



(2) 



where the indices 1,2 refer to the "advanced- retarded" 
decomposition of W, and K = vpnW is the Liouville op- 
erator. This "linearized" action has now the same form 
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as that of a diffusive system, with the diffusion operator 
— DX7 2 replaced by the Liouville operator. Thus, all the 
results derived previously from the linearized action for 
diffusive systems can be directly used for our model, pro- 
vided the eigenvalues and eigenfunctions of the diffusion 
operator are replaced by those of the operator K. 

Since disorder is only present in the close vicinity of 
the boundary, we model it by supplementing the Liou- 
ville operator if by a boundary condition. Generally, the 
boundary condition for an eigenfunction <p(r, n) relates 
its values for outgoing and incoming particles at a point 
on the surface, 



previous motion after the first collision with the bound- 
ary. The system is thus described by two characteristic 
energies, which are the mean level spacing A and the in- 
verse time of flight through the billiard vp/R. This is 
a special type of a billiard, analogous to "hard chaos" 
behavior of genuinely chaotic systems. 

Section VI is devoted to the general situation a < 1. 
A new regime appears for q«1 when the time R/(vpa) 
during which an electron remembers its initial direction 
of motion differs parametrically from the time of flight 
R/vp. In this regime the system is described by three 
distinct energy scales. 



ip(r,n)=\ (Nn')B(n,n')dn' 

\J(Nn')>0 



/ (Nn')B(n,n')ip(r,n')dn', 

J(Nn')>Q 

(Nn) < 0, (3) 



with some kernel B. Here the point r is taken at the 
surface, \r\ = R, and TV is an outward normal to the 
boundary. The form (|j|) ensures that no current flows 
through the boundary of the billiard. 

The scattering kernel B was intensively studied in the 
context of the boundary condition for the distribution 
function (for review, see Ref. ^) and found to be model 
dependent. Realistic models of short-range surface ran- 
domness (like a boundary narrow layer of impurities, or 
a ballistic orifice in a disordered medium) lead to general 
boundary conditions of the type (^), where the kernel B 
is a parameterless function of order one. 

Following Refs. |3l|-|3^, we approximate the above, 
rather complicated, boundary condition by a simpler one, 
where an electron reflects diffusively with probability a 
and specularly with probability 1 — a (0 < a < 1), 



ip(r,n)=na I (Nn') ip(r,n')dn' 

J{Nn')>0 

+ (l-a)<p(r,n"), (Nn) < 0, (4) 

where the vector n" is chosen so that the reflection from 
n" to n is purely specular. Although Eq. (Q) does not 
correspond to any particular known microscopic model 
of disorder, it is commonly believed to provide a good 
qualitative description of surface scattering interpolat- 
ing between purely diffusive (a = 1) and purely specular 
(a = 0) reflection. 

With the exception of SectionlvJ, we consider below 
purely diffusive reflection (Eq. (^) with a — 1). Physi- 
cally, it describes quantum scattering due to short-range 
disorder (correlated at the scale of the order of the wave- 
length 2tt/pp). Alternatively, the diffusive boundary 
condition of this type may result from the purely classi- 
cal scattering off a strongly corrugated surface. As noted 
in the Introduction, atomic-scale disorder has the feature 
that an electron loses the memory about the direction of 



B. Eigenvalues of the Liouville operator for diffusive 
scattering 

As we have mentioned, the level statistics of our bil- 
liard are entirely determined by the eigenvalues of the 
Liouville operator K. Due to the boundary condition (ji|) 
the evolution becomes irreversible, the eigenvalues of K 
have positive real part, and no regularization, like that 
discussed in Refs. II 13 is needed. Below we review prop- 
erties of eigenvalues A of the operator K, 



vpnVip x (r,n) = \ip x (r,n), 
supplemented by the boundary condition 



(5) 




FIG. 1. Natural parameterization of the constant energy 
shell in terms of position r — (r, i?) and direction of motion 
n and its relation to the alternative parameterization by the 
coordinates 8,8',x introduced in the text. 

(p(r,n)=ir (Nn') ip(r,n')dn' , (Nn) < 0. 

J(Nn')>0 



(6) 



The constant energy surface can be parameterized by 
the three real numbers (6,9',x), where the angle 9 (9') 
corresponds to the point where the straight line passing 
through r = (r,d) in the direction n (—n) crosses the 
boundary, and x is the distance from the boundary to r 
along this straight line (see Fig. 111). To cover the whole 
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energy surface the parameters must change in the region 
< 6,6' < 2tt; < x < 2Rsin\(6 - 6')/2\. Eq. (g) then 
expresses the eigenfunction ip at any x through the eigen- 
function at x = 0, i.e. through the function describing 
particles scattered from the boundary. In particular, the 
eigenfunction for the particles arriving to the boundary 
is 



<p x I 6,6',x = 2i?sin 



1 , x = 0) exp sin 

V v F 



(7) 



Now the boundary condition (^J) is used to find a closed 
integral equation for ip\(6,6',x — 0) which is simplified 
by the Ansatz tp x (6,6',x = 0) = (p{6 - 6')exp(iW), I 
having the meaning of the angular momentum, 



<p(0)e 



-iW 



d§ sin0e 2 « sin ^(20), 



(8) 



with the notation £ = RX/vf- Eq. (||) only has solutions 
for the eigenvalues which obey the following equation, 



1 r 

- d6 sin 9 exp [2il6 + 2£ sin i 

2 Jo 



= 0. (9) 



The eigenvalue equation (g) can not be solved ana- 
lytically in a closed form, however the combinations of 
eigenvalues which enter level statistics can be expressed 
through the function J/. Below we list some properties 
of these eigenvalues. 



the relaxation of the corresponding classical system to 
the homogeneous distribution in the phase space. 

The asymptotic form of the solutions of Eq. (||) for large 
|fc| and/or |Z| is given by the saddle-point method: 

J 0.66/ + 0.14 In/ + 0.557rzfc, 0<fc<Z , . 
^ ~ \ (lnfc)/4 + vri(fc + 1/8), 0</<fc' ^ W) 

Note that for k — all eigenvalues are real, while for 
high values of k they lie close to the imaginary axis and 
do not depend on /. Fig. 2 shows a plot of the first 11x11 
eigenvalues of the operator K. 



C. Green's function of the Liouville operator 

The Green's function of the Liouville operator 
T>(ri,ni;r2,ri2) is the time-integrated probability to 
find the particle at the point of the phase space (Yi,ni) 
if it started the motion at (r2,ri2). This function obeys 
the equation 



KV(ri,nr,r 2 ,n2) 
1 



S(ri - r 2 )S(n 1 - n 2 ) - — 



(11) 



where V = ttR 2 is the area of the billiard, and the opera- 
tor K acts on the variables rx,n\. The Green's function 
is related to the correlation of the wavefunctions (Section 
^). It is also used in the alternative derivatio n of th e low- 
energy level correlation function, see Section III B . 
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FIG. 2. The first 11 X 11 (0 < k, I < 11) eigenvalues of the 
Liouville operator K in units of vf/R, as given by roots of 
Eq. §. 

For each value of / = 0, ±1,±2, . . . Eq.@ has a set 
of solutions £,ik with ^ = £,-i.k = Cf-fe: which can 
be labeled with k — 0, ±1,±2, ... (even I) or k = 
±1/2, ±3/2, . . . (odd I). Thus, the eigenvalues form a 
two-parameter set. For / = k = we have £oo — 0; cor- 
responding to the zero mode ip(r,n) — const. All other 
eigenvalues have positive real part Re^fc > and govern 



1. Green's function integrated over momenta 

In this subsection we calculate the Green's function 
integrated over momenta: 



ns(ri,r 2 ) = J dnidn 2 X>(ri,ni;r 2 ,n 2 ), 



(12) 



which describes the probability to find an electron at r\ 
after it has been found at r 2 , and is important for the 
wavefunction correlation. Integrating Eq. ([ll]) over <in 2 , 
we obtain 



vf n-i 



dh 
dri 



1 

TIP 



S(r 1 -r 2 ) - — 



(13) 



where /i(ri,ni;r 2 ) = J G?n 2 2?. To solve Eq. (13), we 
use the same strategy as with Eq. (|^) and replace the 
coordinates (rx,ni) by the variables 6i,6[,xi, 



h{6 l ,6' l ,x l ;r2) = h{6 1 ,6' x ,x l = 0;r2) 



2 r Xi 

+ — / dx'l 

VF Jq 



5(r(6 



i,o 1 ,x 1/ 



r 2 



1 

V 



(14) 



The function h(6i,6' 1 ,xi = 0;r 2 ) which describes the 
particles leaving the boundary at the point (i?, 9\ ) , does 
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not depend on 8[ . Using the boundary condition for the 
operator K, we obtain the equation for this function (be- 
low the irrelevant arguments are dropped and the func- 
tion is denoted by h(8i;r 2 )), 



h(6i;r 2 ) 



d8" sin — h{9 1 



271 



2Rsin(\e 1 -8"\/2) 



. d9" 

2 Pf Jo jo 



dx'l sin — 
1 2 



1 

V 



(15) 



Calculating the second term in the r.h.s., we obtain 

h(8 1 -r 2 ) = - f dS"sm^h{ex+e";r a ) (16) 
4 Jo z 



1 

PF 



F(r 2 ,9i - 0a) - 77 / r 2 dr 2 dti 2 F(r 2 , t - tf 2 ) 



Here we have introduced the function 

Fir 8) = R - r COS 9 
l7 *' ' R 2 + r 2 - 2i?rcos6»' 

and used polar coordinates (rk,"&k) f° r the point r&, 
fc=l,2. 

Eq. (|l6| ) is solved by expanding the function h{6\, r 2 ) 
in a Fourier series; in this way we find all the Fourier 
components hk except for fc = 0. The component ho is 
not determined by Eq. ([n3|) and must be fixed by the 
conditions that IIb(t"i, v 2 ) is symmetric and its integral 
over dr\ equals zero (the conservation of the number of 
particles). Restoring subsequently the function h in the 
bulk of the billiard from h(8\;r 2 ) by means of Eq. ( |l4| ) 
and integrating it over dn\, we obtain the Green's func- 
tion of the Liouville operator integrated over momenta, 



n B (ri,r 2 ) = /o(ri,r 2 ) + /i(ri,r 2 ), 

J- J* - 

ttpf [ |n - r 2 \ V 

1 /" 1 1 

V 



fo(ri,r 2 ) = 
df 2 
fi(rx,r 2 ) = 



ri - r 2 \ 
1 



V 2 
,4k 2 



d?i- 
d?idr 2 -r^ 



n 
i 



r 2 \ 



k—l 



1 ( r x r 2 
4fc 2 V R 2 



ri - r 2 \ 



COS fc(#l 



(17) 
i? 2 ). 



It can be actually traced in the course of the calculation 
that the term fo in Eq. ( |l7| ) describes the propagation of 
an electron from r 2 to ri along the straight line (no col- 
lisions with the surface). Likewise, the term /i describes 
the processes which involve at least one collision; further- 
more, in the factor 4fc 2 — 1 in /j the term 4fc 2 originates 
from the trajectories with one collision, whereas —1 is 
related to the double and multiple collisions. 



2. Full Green's function 

Analogously to the previous subsection the value of the 
function T> in the bulk of the sample is expressed through 
its value on the boundary. Instead of Eq. ( |l4[ ) we now 
have 



V(9i,6' 1 ,xi;r 2 ,n 2 ) 
~ ' dx'l 



PF Jo 



X>(0i,0i,a:i = 0;r 2 ,n 2 ) 
Sir^e'^x'l) - r 2 )5 {m - n 2 ) - 



V 



(18) 



Due to the boundary condition, the function 
T>{6\\ r 2 , n 2 ) = V(0i,6' 1 ,xi = 0; r 2 , n 2 ) does not depend 
on and obeys the equation 

X>(0i;r 2 ,n 2 ) = i f d6"sin -^V(6 1 +6"-r 2 ,n 2 ) 
4 Jo * 



1 

2p~F~ 



d6" 



2_Rsin(|6»i-e"|/2) 



dx'l sin ■ 



S(r(0",e 1 ,x'{)-r 2 )S(n(9",0 1 



n 2 ) 



(19) 



After lengthy calculations, we find that the second term 
in the r.h.s. equals (2tt/pfR)[5(8 1 — 8' 2 ) — 1/2tt], where 
8' 2 is the polar angle corresponding to the point at the 
boundary to which the vector n 2 points from r 2 . For 
definiteness, we restrict all the angles 6-y, 8[, 8 2 , 8' 2 to the 
interval [0, 2ir]. Continuing in the same way as in the pre- 
vious section, we find the final expression for the Green's 
function of the Liouville operator, 



T>(ri,m;r 2 ,n 2 ) =V +Vi 
1 f , f ri - r 2 



(20) 



-Do = 



Trp F 
1 

V 



77 / d ?i S 



ri - r 2 \ 
fi - r 2 



- n -\ 



- dfoS 



T>i 



1 

V 
1 

V 2 
2tt 



\ri 



r 2\ 

f 2 



Til 



ri 

df\df 2 S 



Rpi 



ral 


m 


fi - f 2 


fi - f 2 \ 


1 " 




~ 2^ 


i 



5(n 1 - n 2 ) 
\ri - r 2 \ 
5(ni - n 2 ) 
r 2 \ 
- n 2 ) 



\ri- 
8(m 



ri - r 2 \ 
5(ni 



Til 



n 2 



ri - r 2 \ 



T>o = 



24Rp 



5(6! - 8' 2 ) 

[3(0i - 8' 2 ) 2 - 6tt|0i - 



2tt 2 



It is straightforward to check that the integration of Eq. 
( pp| ) over drt\dn 2 gives Eq. (0). The Green's func- 
tion (20) possesses obvious properties T>(ri, n±; r 2 , n 2 ) = 
T>(r 2 ,-n 2 ;rx,-nx) and J dr 1 dn 1 V(r 1 ,n 1 ;r 2 ,n 2 ) 
| dr 2 dn 2 V(ri,n 1 ;r 2 ,Ti2) = 0. 

Actually, T>o is a solution of Eq. (|ll|) in the infinite 
space and thus represents propagation processes which 
do not involve any collisions with the boundary. This 
kind of propagation is only possible if ri\ coincides with 
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ri2 and both are directed along r% — r%- Furthermore, 
T>\ is responsible for the propagation from r2 to r\ with 
only one intermediate collision, which necessarily requires 
9 1 = 9' 2 . Finally, the term T>2 describes propagation 
which involves two or more intermediate collisions. 



III. SPECTRAL STATISTICS 

In this Section we discuss the level correlation function 
and level number variance of the circular billiard with the 
boundary condition (||). 

As we have mentioned in the Introduction, the results 
for the level statistics follow readily from general formu- 
lae derived for the diffusive conductorso where the eigen- 
values of the diffusion operator are replaced by those of 
the Liouville operator. 



A. Low frequencies 

We define the level correlation function in a standard 
way, 

R 2 {u) = {VA) 2 {v{e + u)v(e)) - 1, 

where v(e) is the (fluctuating) density of states, A = 
{Vv)~ l is the mean level spacing. 

In the range of relatively low frequencies (which in our 
case means u> <§; vp/R, seeubelow) the function i?2(w) 
quite generally has the formEfl (s — uj/A) 



Si{w) = -(R/vf? 



VF 

RA 



sin 7rs 
p F R 



A . 2 

„ „ sin its, 



(21) 



The first two terms are universal and actually correspond 
to the random matrix theory result. They are associ- 
ated with the contribution of_the mode with zero eigen- 
value of the Liouville operator!! The last term represents 
the non-universal correction; the information about the 
operator K enters through the dimensionless constant 
A = Yli ^hl i wfiere the prime indicates that the eigen- 
value £oo = is excluded. The value of A, as well as 
the high-frequency behavior of i?2 (sV-feee below) , can be 
extracted from the spectral functionESl 

S{uj) = Y j Si{lo); Si(Lu)=Y,^ki-itor 2 ■ (22) 



According to the Cauchy theorem, Si can be represented 
as an integral in the complex plane, 



Si(uj) 



R\ 1 



1 



v F ) 27Ti J c (z - iujR/ v F ) 2 Ji (z) 



dz, (23) 



where the contour C encloses all zeroes of the function 
Ji(z). Evaluating the residue at z = iuR/vp, we find 



dz 2 



In Ji(z). (24) 



■jR/vp 



Considering the limit uo — > and subtracting the contri- 
bution of A o = 0, after lengthy but elementary algebra 
we obtain 

= -19/27 - 175tt 2 /1152 + 64/(9tt 2 ) w -1.48. (25) 

In contrast to the diffusive caseEi, this constant is neg- 
ative: the level repulsion is enhanced with respect to re- 
sult for RMT. We recollect that in the diffusive case the 
level repulsion is suppressed as compared to RMT, and 
this suppression has a plausible physical explanation. In- 
deed, the non-universal correction in diffusive case is pro- 
portional to g^ 2 , where = E c /A 3> 1 is the dimen- 
sionless conductance, and E c is the Thouless energy. In 
disordered conductors, the parameter gd is responsible 
for the metal-insulator transition/crossover. This non- 
universal terrn-thus reflects a tendency to localization in 
level statistics^ and would become of the order of unity 
for gd ~ 1, when the system approaches the insulating 
regime (with uncorrelated levels, Rzfs) = S(s)). 

In contrast, the last term of Eq. (21 ) describes different 
physics. For ballistic systems, the limit gt ~ 1 means not 
the insulating regime, but rather a quantum-mechanical 
system far from the semi-classical regime. Energy levels 
in these systems may be strongly correlated, even more 
strongly than in RMT, hence the tendency to enhance- 
ment of the level repulsion observable from Eq. ((2^) . We 
will see that this tendency is even more pronounced in 
the high-energy behavior of the level correlation function. 

The last point we make is the following. Eq.(fj"l|) is 
valid as long as the correction is small compared to the 
RMT result, i.e. provided ui is below the inverse time of 
flight, vf/R- Thus, the inverse time of flight serves as a 
Thouless energy for this ballistic billiard, with gi, play- 
ing the role of the dimensionless conductance. Note that 
gb is related to the number N of levels below the Fermi 
energy as gt, — N 1 / 2 . 



B. Low frequencies — Alternative derivation 

The constant A can also be found in another way, using 
the Green's function of the Liouville operator. Indeed, 
the quantity S(0) can be written in the following form, 

S(0) — (irv) 2 J dridr2dnidn2T>(ri,ni;r2,ri2) 

xV(r2,n 2 ;r 1 ,n 1 ), (26) 



where the function V is given by Eq. (poj). Eq. (|26|) is 
directly verified by expanding the Green's function in the 
eigenf unctions of the operator K. Now we calculate the 



G 



integral in Eq. ( p6| ) directly. It is instructive to split the 
result into three terms, 



5(0) = S A + S B + S c , 
Sb = J (p V! + V{D Q + VtV x ) 

8 7T 2 ' 

S c = J (V T>2 + V 2 V + V{D 2 + V 2 V X + V 2 V 2 ) 



R_ 

Vp 



(27) 



which reproduces Eq. (p5|). In our notation, 5^ is 
given by the processes which involve no collisions with 
the boundary, Sb is a contribution of trajectories with 
only one collision, and Sc takes into account all other 
processes. We see that all these contributions are of the 
same order and can not be disregarded. 



C. High frequencies 

In the range u> ^> A the level correlation function can 
be decomposed into the smooth (Altshuler-Shklovskii) 
parto, 



R^ S (lo) = {A 2 /2n 2 )ReS(co), 



(28) 



and the part-J?2 SC which oscillates on the scale of the 
level spacingp3. We consider below the high-frequency 
regime w ^> vp/R. 



1. Smooth part 

For high frequencies, the integral in the definition (j|) 
of the functions J; is small, and thus J; is close to —1. 
Expanding the logarithm in Eq. ( p4| ) up to second order 
in Ji + 1 and using the Poisson formula, 



^e 2ae =TT V 6(9 -7m), 



(29) 



to sum over I, we obtain for lj 3> vp/R 



Vp 



Thus, the smooth part of the level correlation function 
is an oscillating function of frequency, with the period 
of itvf/2R and a slowly decaying amplitude. To clarify 
the connection with the periodic orbit theory, we note 



that any closed sequence of chords joined at vertices is a 
"periodic orbit" in our model. The above period corre- 
sponds to the "periodic orbit" which traverses twice the 
diameter of the billiard (the longest "periodic orbit" with 
two vertices). The amplitude of the oscillations is pro- 
portional to g b 2 ■ This is very different from the diffusive 
regime, where the smooth part of the level correlation 
function does not exhibit any oscillations. Furthermore, 
in 2D case the AS contribution vanishes, and theJeading 
behavior is provided by weak localization effects^!!. 



2. Oscillating part 

The oscillating part of the level correlation function 
i?2 SC ( s ) f° r frequencies to A is given byE3 

^2 SC ( S ) = (l/27r 2 )cos(27rs)L>(s), (31) 

where D(s) is the spectral determinant, 

D(s) = S - 2 J] (l^sA/A^r^^sA/A,,)- 1 . 

fci^(OO) 



(32) 



Since A- 2 d 2 \nD(s)/ds 2 = -2ReS(aA, the spectral 
determinant can be recovered from Eq. ([24]), 



D 



<•> - (I 



6 p ci+c 2 s 



n 



i 



Ji(is/g b )Ji(-is/g b ) ' 



(33) 



where c± and c 2 are arbitrary constants which are fixed 
by the requirement that Eq.(|3l]) in the range A«w< 
vp/R reproduces the low- frequency behavior (|2l|). Ex- 
panding InD in s/g b <C 1 and comparing it with the 
low-frequency expression 



lnD(s) 



1 



2 ' 



which follows from Eq. (f2l|) , we obtain c\ = c 2 =0. Fi- 
nally, for s > 1 (uj > vp/R) Eq. (||) yields for the 
oscillating part of the level correlation function, 



128 5b 2 



2ttlo 



(34) 



It is remarkable that the amplitude of the oscillating part 
does not depend on frequency. This is in contrast with 
the diffusive case, where in the AS regime (u> above the 
Thouless energy) the oscillating part R 2 sc (lo) is exponen- 
tially smallE3. Such a behavior indicates a strong rigidity 
of a level system and is reminiscent of ID harmonic os- 
cillator. 
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D. Spectral form- factor 

To illustrate the nature of the oscillating terms in 
the asymptotes (|o|), (p^) we compute numerically the 
form factor K(t) = J due luJT R%{uS) of the two-level cor- 
relation function for gb = 5. The result is shown on 
Fig. [| The non-zero limit at to — ► oo of the oscillat- 
ing part i?2 SC (w) gives rise to a 6- function contribution 
of the form TT 5 /^28g 2 )5(T — 2ir/A) (shown as a vertical 
spike on Fig. [|)tia. The "periodic orbit" which traverses 
the billiard twice along the diameter and has the pe- 
riod Ti ~ AR/vp manifests itself as a set of power-law 
singularities in K(t) at r = T 2 and r = 2ir/A ± T 2 . 
Near to these singularities the form factor diverges like 
1/VT 2 - t and l/y/T 2 ± (2tt/A- r), respectively. The 
peak to the right of the first singularity is the contribu- 
tion from the "equilateral triangle" periodic orbit with 
the period T 3 = 3y/3R/v F . The same orbit causes small 
peaks at r = 2n/A±T 3 . 



<5(t-2tt/A) 




FIG. 3. Form factor for the two point correlation func- 
tion K(t) = J du>e*" r i?2 {&) for gb = 5. The 5-function at 
t = 27r/A is shown as a vertical spike. It arises because of 
the non-decaying oscillations (^). The 1 / s/ T2 — t singular- 
ity at T2 = 4R/vf is related to the "periodic orbit" which 
has period T2 and traverses twice the diameter of the billiard. 
The same orbit gives rise to singularities at r = 2tt/ A ± T2 



= / (s - \s\) R 2 (s)ds 

J—s 



(35) 



Direct calculation gives for s <C gb (E <C vf/R) 

As 2 

tt 2 E 2 (s) = 1 + 7 + ln(27rs) + — - 

^9b 

and for s > <?& (15 » vf/R) 
7T S 2 (s) = 1 + 7 + In — 5- 



(36) 





Us 












\9b 





(37) 



Here 7 w 0.577 is Euler's constant, and ^4 is define d by 
Eq. (P5[) . The first three terms at the r.h.s. of Eq.(|3q) 
represent the RMT contribution, which is logarithmic in 
energyEl It is remarkable that the actual level number 
variance is always smaller than that given by RMT. 



E 2 (s) 



JO) 




s/g b 

FIG. 4. Level number variance £2 (E) as a function of en- 
ergy; s = E I A. Curve 1 shows the RMT result, while curves 
2 and 3 correspond to asymptotic regimes of low Q3q) and 



high (|3 7|) frequencies. The saturation value E 2 
Eq. (jm 



(()) 



is given by 



E. Level number variance 

In practice, the description of level statistics by the 
means of the level correlation function is not always con- 
venient. Indeed, for low frequencies uj <C vf/R non- 
universal effects are manifested only as small corrections 
to RMT. For high frequencies, the entire behavior is non- 
universal, but in this range the level correlation function 
is proportional to g^ 2 and small. Therefore, in order to 
study non-universal effects in level statistics, experimen- 
tally or by means of computer simulations, it is instruc- 
tive to consider a quantity more sensitive to these effects. 

A well-known way to emphasize the high-energy be- 
havior of the level correlation function is to study the 
variance of the number of levels in an energy interval of 
width E = sA, 



As seen from Fig. |], the two asymptotes ( p6| ) and 
d37| ) perfectly match in the intermediate regime, s ~ gb 
(E ~ Vf/R)- Taken together, they provide a complete 
description of £2(5)- According to Eq.(^), the level 
number variance saturates at the value 



-.(0) 



2 (l+7 + ln(I6<7 6 /^ 2 )) > I- 



(38) 



This saturation again confirms the conclusion that we 
have already made from analyzing the low-energy behav- 
ior of the level correlation function — the system of levels 
of our billiard is quite rigid, more rigid than in RMT. We 
also remark that a saturation of E 2 (s), as well as its os- 
cillations on the scale set by the shortest periodic or 
is predicted by Berry for a generic chaotic billiardE 
This behavior is also in agreement with the results for 
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£2(5) found numerically for a tight-binding mcyjd with 
moderately strong disorder on boundary sitesEJ£3tl 

At the same time, we note that the behavior of £2(3) 
is quite different from that in othpc reference systems. 
Indeed, for both diffusive systeHtsEijand systems with 
bulk disorder in ballistic regimetlnla the level number 
variance is greater than RMT. The energy at which £2 
is expected to saturate depends on the type of disorder. 
For short-range impurities (white noise random poten- 
tial) arbitrarily short periodic orbits exist, and thus no 
saturation up to u ~ Ep is expected. For a diffusive sys- 
tem and smooth random potential the diffusive dynam- 
ics lead to a linear increase of the level number variance, 
£2(5) ~ s/gd, for s 3> Qd ~ vvpltr, while the shortest or- 
bits have lengths of the order of the transport mean free 
path lf r , causing the saturation of £2(5) at a paramet- 
rically larger value s ~ vp/(Z tr A). Similarly, in-ikough 
billiards" (slightly distorted integrable billiardaid'EJ) the 
level number variance is higher than in RMT and does 
not saturate until a value s ~ vp/RA parametrically ex- 
ceeding the effective Thouless energy, since the system is 
diffusive in the angular momentum space. 



IV. PARAMETRIC LEVEL STATISTICS 



A. Introduction 

In this Section, we study the parametric statistics of 
energy levels of our system. Specifically, we assume that 
the billiard is placed in a magnetic field, which plays the 
role of an external parameter. 

Alreadya considerable amount is known on the subject 
(see Ref. |42| for review). In this paper we are interested 
in the parametric level correlation function 

R^{uj,B) = -\ (39) 
+ (VA) 2 (u(E + cu/2, B + B/2)v{E - w/2, B - B/2)) , 

where the mean magnetic field B is introduced in or- 
der to break time-reversal symmetry. The correlation 
function ( |39| ) has beery-oreviously investigated in .the di- 
agrammatic expansions, the cr-model approachc.sE3, the 
random matrix-jLheory (by means of Dyson's brown 
motion model)c3o, and by semi-classical methodscJ 
The most remarkable observation of these works is that 
when the perturbation is weak and the frequency u> is 
low, the parametric correlation function is universal, i.e 
it does not depend on the type of the perturbation, pro- 
vided it is properly rescaled, nor on any details of the 
system. 

Below we study non-universal behavior of the paramet- 
ric level correlation function (|39|). We use the cr-model 
formulation of parametric statistics developed in Ref. 
to derive an analog of Altshuler-Shklovskii formula for 
this function and evaluate it explicitly for our model. 



B. Eigenvalues of the Liouville operator in magnetic 

field 



The correlation function ( p9[ ) may be obtained from 
the supersymmetric cr-model with the effective action 
(|l|) provided the operator V is replaced by the "gauge- 
invariant" combination V — (ie/2c)AA, where A is the 
vector potential which corresponds to the field B: B = 
V x A. Repeating the steps leading to Eq. (g) we find 
that in the effective action the Liouville operator K is re- 
placed by its gauge-invariant form in the magnetic field, 



Kg, = v F n ( V A 



(40) 



Prior to the investigation of statistical properties of the 
energy levels, we find the eigenvalues of the operator K$. 
Choosing the symmetric gauge 4 = B x r/2 we obtain, 
instead of Eq. (||) , the equation 

Jf (0 = -1 (41) 

1 r 

+ - d0 sin 9 exp [2iW + 2£ sin 9 + i<\> sin 29] = 0, 

2 Jo 

which determines the dependence of the eigenvalue A; = 
^vf/R on the magnetic field. We have introduced the 
dimensionless parameter </> = $/$0) where $ = ttBR 2 is 
the magnetic flux through the billiard, and $0 — 2?:hc/e 
is the flux quantum. 

Equation ( |i~L| ) defines a two-parameter set of eigenval- 
ues £ki- It stays invariant under the simultaneous trans- 
formations I —t — I, 4> — 4>, which means that 
if X(B) is an eigenvalue of the operator K§ then A*(— B) 
is also an eigenvalue. 

For sufficiently low frequencies and magnetic fields the 
parametric correlation function (^) is dominated by the 
zero mode I = k = (see below) . To find the evolution of 
the corresponding eigenvalue £00 , we put I = and take 
a variation of Eq. ( ffl| ) with respect to </> <C 1, 



d6 sin6» 



2£oo sin 9 + i<j> sin 29 - y sin 2 29 



0. 



which gives £00 = 8(j) 2 / (15n) . 



C. Parametric level correlation function 

The parametric correlation function i?$ ( |39| ) is ex- 
pressed in terms of the retarded G R and advanced G A 
Green's functions in the following way, 

R*(w,B) = ±[T(u,B)+T(-u,-B)], (42) 



T(cu,B) 



2{ttv)' 



(Tr G R (E + uj/2, B + B/2) 



TtG a (E-uj/2,B-B/2)) c 







where (AB) C = (AB) — (A)(B) denotes the irreducible 
part (cumulant). 

For low frequencies and magnetic fields (the precise 
condition is specified below) the function i?$ is non- 
perturbative and can be found using the random ma- 
trix theory. This regime was previously investigated in 
Ref. 



44 



Here we focus instead on the case of higher 
frequencies and/or fields, when the smooth part of the 
parametric correlation function R^ s (cu,B) is correctly 
described by perturbation theory. From Eq. ( [42] ) we 
obtain 



r£ s (u,b) 

A 2 x 

kl 



4tt 2 

A 2 
2tt 2 



1 



ReV- 

^ \— i. 



-iuj + Xu(B)} 2 
1 



[ioj + \ kl (-B)f 



k-l 



iw + A w (B)] 5 



To identify relevant parameters, we consider R§, s at 
zero frequency. For low fields, the sum in Eq. (|43| ) is 
dominated by the lowest eigenmode, which as we have 
seen is quadratic in field, Aoo ~ vf4> 2 /R- The per- 
turbation theory does not apply for Aoo IS A, which 

corresponds to the fields 4> g b 1 ^ 2 , where as before 
.96 = vp/(RA). Generally, a non-perturbative calcula- 

1/2 

tion is needed when ui ^$ A and <fi < g b (region 1 in 
Fig. |^) . Everywhere outside this regime, the perturbative 
expression (|43|) applies. 



<fi ~ 1. At <f> ~ 1 other modes become important, and, 
in addition, the zero-mode eigenvalue can not be taken 
quadratic in field any more. Thus, the point <p ~ 1 plays 
for parametric correlations the same role as u) ~ vp/R 
for the usual (non-parametric) level correlation function. 
Generally, for <f> -C 1 and u> <C vf/R we are in the 
regime when everything is determined by the zero-mode 
approximation. For higher fields <fi 3> 1 or frequencies 
w Vf/R, the system crosses over to the Altshuler- 
Shklovskii regime (region 3 in Fig. 5), when all the modes 
are important. 

To end our qualitative discussion, we determine the 
strength of the magnetic field needed to strongly affect 
the classical dynamics. The cyclotron radius of an elec- 
tron trajectory in the magnetic field is r c = mvpcj (eB). 
The magnetic field strongly affects the dynamics pro- 
(43) vided r c < R, which gives <fi > g b . Our theory is thus 
valid for <f> <C g b , which still leaves a large window for the 
regime <f> 1 . 

1. Perturbative zero-mode regime (Region 2 in Fig. 5) 

Using the result for the zero-mode eigenvalue, we find 
fromEq. (|§) 



-1/2 




A v F /R to 

FIG. 5. Different parameter regions (in the frequency 
magnetic field plane) for the parametric level statistics: 1 
non-perturbative (RMT) region, 2 - zero-mode region, 3 
Altshuler-Shklovskii region 



R 



AS, 



225 64 5 2 4 - 225tt 2 s 2 
2 (64 5 2 4 + 225tt 2 s 2 ) 2 ' 



(44) 



where, as previously, s = ui/A. The result (03) applies 

— — 1/2 

when both <p <C 1 and s <C g b , but cither <fi g b or 
s ^> 1. Without the magnetic field, </> = and we repro- 
duce R$ s (s,0) = — (2tt 2 s 2 )~ 1 , which is a smooth version 
of the RMT result (Ell). For s = 0, we obtain 



i?£ S (0, 



225 



2. Altshuler-Shklovsktt regime (Region 3 in Fig. 5) 



For <|i > 1 or s > 3(, we generalize Eq. (24) to the case 
of parametric statistics, 



Ri s (s, 



■ \ „ ReT^ 



dz 2 



s/9b 



lnjf(z), (45) 



For higher fields, the zero mode still dominates until + 1 up to the second order, we obtain 



and use the fact that J[ 



-1. Expanding In J l f(z) in 



R£ s (s,cf)) = -— -L^Re j~ 2 ^ d0 sin0exp(2iZ0 + 2*sin0 + #siii.20) 
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i d 2 r 

--t-j / e$i e$2 sin 0i sin 2 exp [2il(6x + 2 ) + 2z(sin0i + sin0 2 ) +#(sin20 a + sin20 2 ) 
8 dz Jo 



(46) 



z=is/g b 



Now we employ the Poisson formula (|29|). The first 
integral in Eq. (46) vanishes identically, and in the sec- 
ond one we have 9\ + 02 = 7J". Under this condition, the 
integral does not depend on the magnetic field. We thus 
have 

1 f n 

i?£ s (s,0) = —j Re / d0 sin 4 0exp(4issin0/ 5b ) (47) 

f (8 5h 2 r 1 (3-20 S 2 /^ 2 ): 
~ \ (2vrs.g3)- 1 /2 CO s(4s/ 5fc - tt/4), s > 5b ' 

The asymptotes in the regime s 3> g& are identical to Eq. 

The result that the parametric level correlation func- 
tion does not decay with magnetic field is very different 
from the diffusive caseEJ and is actually quite surprising. 
It means that the density of states at the same energy re- 
tains the memory even at very large fields, after the levels 
underwent many avoided crossings. The field dependence 
of the correlator appears if we keep higher orders of the 
expansion of the logarithm in Eq. ([45|), but it is only 
manifested as small corrections to the result ([47|). 



D. Parametric level number variance 

The high-frequency behavior of the non-parametric 
level correlation function i?2(s) is best seen in the level 
number variance £2(3), see Section [II. Similarly, it is in- 
structive to introduce the parametric level number vari- 
ance (PLNV)B, 



U{E, B) = ( [N(E, B) - N(E, B + B)\ 



(48) 



where N(E, B) is the full number of levels in the energy 
strip E in the magnetic field B. In zero field PLNV van- 
ishes (in low fields it is generally linear in field, see Ref. 
g). If the values of N(E,B) and N(E,B + B) are uncor- 
rected in high fields (which is not the case in our situa- 
tion, since the level correlation function does not decay in 
high fields), PLNV saturates, U(E, B -> 00) -> 2E 2 (E). 

PLNV can be easily expressed through the level corre- 
lation function (analogously to Eq. (35)), 



U(s, <f>) = 2 (s- \S\) [R 2 (S) - R*(s, (/>)] dS, (49) 

J — s 

where s = E/A. Below we concentrate on the perturba- 
tive regime <f> 3> g b 1 ^ 2 , when 



(s — |s|)i?$(s, 4>)ds 

J<l>l 



-^ReJ] [in. V(,.s /,,,,) In./; (in 
I 



(50) 



— 1/2 

A direct calculation gives for g b <C (j) <C 1 

U(s,<t>) 



2 N 225 s 2 

— (In27rs + I + 7) 



u(s,4>) 

U(s,<t>) 



7T 



7T" 

2 



128 5 2 ^ 4 ' 



1< s < g h tf 

In h 1 + 7 

15 

1 32gb ^ 1 1 

In h 1 + 7 

15 ' 



(51) 



g b <t> < s < g b , 
-A', s><? 6 , 



where 7 is again Euler's constant, and 
n 2 l 4/ 2 - 1 



0.115. 



In higher fields 4> ^> 1 we have 



C/(s,0) = ^(ln2vr.s + 1 + 7) + 



9l W 2 



C/( s ,0) 



1< s < g 6 , 

1 16 # b , 1 , ' 
m — 5- + I + 7 



(52) 



1 



where the constant A is given by Eq. ( p5|) . A comparison 
of Eq. @ with Eqs. @ and @ shows that in the 
limit ^ — > 00 PLNV differs from 2Ea(s). This difference 
originates from the fact that the correlation function i?$ 
does not vanish at large 4> (see Eq. (E7l)). 



V. CORRELATIONS OF EIGENFUNCTIONS 

According to Berry's conjecture^, a wave function in 
a 2D chaotic system shows Gaussian fluctuations with 
the correlation function V(ip*(r)i/j(r')) = Jo(ppr) (see 
Ref. |5(] for a recent generalization of this conjecture). 
The supersymmetry method allows one to derive this re- 
sult (which is equivalent to the zero-mode approxima- 
tion for thx^ (jLjgodel) and to calculate system-specific 
corrections^TfJO . 

Let us consider the two-point correlation function 

C(r 1 ,r 2 ) = A^|^(r 1 )^(r- 2 )|5(e- ea )\ . (53) 

In the standard diffusive <7-modelE], the correlation func- 
tion (^3|) can be written in the form of the integral over 
the c-model field Q(r), 
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riA 

C(r 1 ,r 2 ) = lim (Gn(ri, ri)G 22 (r 2 , r 2 ) 

+ Gi 2 (ri,T- 2 )G ! 2i(T-2,ri)) S [Q] , 



(54) 



where r/ = —iuj/2 > Ois the level broadening, G is 
the Green's function in the field Q, and the subscripts 
1,2 refer to the advanced- retarded decomposition, the 
boson-boson components being always implied (we drop 
the corresponding indices). Eq. ( |54] ) yields the following 
resultEj, 



V 2 C(r l7 r 2 ) ~ 1 + n D (n, r 2 ) + fc,(n - r 2 ) 
'r 1 +r 2 r!+r 2 " 



-n 



D 



k q (ri - r 2 ), 



(55) 



propagator, k q (r) 



where IId is the diffusion 

Jq (pf r)e~ T l l 1 and Z is the mean free path. In the 
framework of the ballistic cr-model the diffusive propa- 
gator Il£)(ri,r 2 ) is replaced by its ballistic counterpart 
IlB(ri, r 2 ) given in our model by Eq. (|T7j). At this point 
we seem to encounter a problem. Indeed, at short dis- 
tances r =\r\ — r 2 \ the classical propagator ils(ri,r 2 ) 
is dominated by the contribution of the direct path, 



n s (ri,r 2 ) ~ / (ri,r 2 ) 



1 



npF \ri - r 2 \ 



(56) 



This contribution, which becomes of order unity at r ~ 
Xf, would imply strong deviations from the universal 
Gaussian statistics. However, it is not difficult to re- 
alize that this contribution is nothing else but a clas- 
sical "copy" of the term k q (ri — r 2 ). Therefore, we 
encounter a problem of the double counting: one and 
the same contribution is taken into account twice, classi- 
cally and quantum-mechanically. In the case of the dif- 
fusive cr-model this problem does not appear because of 
the scale separation: the classical propagator II d is re- 
stricted to low momenta q < l^ 1 , while the short-scale 
(r < I) physics corresponding to high momenta (q > 
is treated quantum-mechanically. 



The situation is different in ballistic case. The semi- 
classical description extends now to all momenta q < pf- 
Therefore, there is no separation in momentum space be- 
tween the slow modes (treated within the semi-classical 
approximation) and the fast ones (treated exactly, i.e. 
quantum-mechanically), and a careful treatment is re- 
quired in order to avoid the double counting. For this 
purpose, we find it instructive to consider a problem 
with an additional smooth random potential U(r) in the 
bulk, characterized by a correlation function V(i — r') = 
(U(r)U(r')) with a correlation length d^> \p and induc- 
ing a small-angle scattering. Specifically, we will assume 
that the corresponding transport scattering rate r^ 1 is 
negligible, i.e. the transport mean free path l tr = vpTtr is 
large, k r S> R, so that the bulk scattering has essentially 
no effect on the classical dynamics. On the other hand, 
the single-particle mean free path I — vft (corresponding 
to the total relaxation rate r" 1 ) will be assumed to sat- 
isfy Xf "C d -C I -C R and will play a role similar to that 
of ltr for the diffusive cr-model, separating the regions of 
classical and quantum treatment. The condition <i <C i 
ensures that the scattering on this random potential is of 
quantum-mechanical (rather than classical) nature and 
is correctly treated within the Born approximation. 

We first ignore the boundary scattering and will return 
to it in the end of the calculation. Averaging over the re- 
alizations of the smooth random potential U(r), one can 
derive the cr-model following Ref. ^5|. (This derivation is 
outlined in Appendix). 

To calculate deviations gf-thejjeigenfunction statistics 
from universality, we writeQoHo (See Eq. (A4)) 



T(r,n) = T (l-W(r,n)/2) 



(57) 



and then integrate out perturbatively non-z ero modes de- 
scribed by W(r,n). The part of action (AS) which is 
quadratic in W has the following form in the momentum 
space, 



.Fo = T^Str J dn 1 dn 2 (dq)W 2 i(-q,ni)Wi 2 (q,n 2 ) 



< I— g J {dp)w(n 1 ,n)w{n,n 2 )G R {p + )G A {jP-) +irvw(n ll n 2 ) 
^Str J dnidn 2 (dq)W 2 i(-q,ni)Wi 2 (q,n 2 )[-A(q]ni,n 2 ) + B(ni,n 2 )} , 



(58) 



where p± — p ± q/2. In the last line of ( |58|) we intro- _ 2V(q-n 2 raOStrP - — — R - ^ 

duced the definitions A and B for the two terms of the " 2 2 

quadratic form. Eq. (pm induces the contraction rules . \ c , n l + A 1 — A 

t . , , M ytLJ , w , , „ f +2D(—q; m, n 2 )StrP — - — R — - — 

tor integrals over non-zero modes W (cj. Ref. Kq), 2 2 



Str[W(-Q, ni)P]Str[W(q, n 2 )R] 



StrW(-q,ni)PW(q,n 2 )R 
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= 2V(q;n 2 ,n 1 )StxP^A StrR 1 A 



-2£>(-g;ni,n 2 )StrP 



2 

1- A 



StrR 



2 

1+A 



(59) 



where P and R are arbitrary matrices and the propagator 
T> is given by the series 



V = {-A + B)- 1 = B- 1 + B- x AB~ l 



(60) 



We are now ready to evaluate the ballistic counterpart 
of Eq. (Q) (with the cr-model field being Q(r,n)). The 
Green's function G there is given by 



C [E-Hq- - I dn'Q(r, n')w(n, ri) ) (61) 



Expanding the term (G11G22) in Eq. ( p4[ ) up to 
quadratic order in W and switching to the momentum 
space, we get 

J d(n - r 2 )<Gii(ri,ri)G a2 (r2, r 2 )>e*< Pl - pa > 

= -i/(dp l)rf n>(n l! n' 1 )(dp 2 ) ( in>(n 2 ,n 2 ) 

x([G r o(Pi_)<3i(-9,ni)G (Pi+)]n 
x[Go(p 2+ )Qi(q,n' 2 )G (p 2 _)} 22 ) , (62) 



where G = Htd.Gr ■ Q + ReG R , Q = T AT , and 
Qi(q,n) = T AW(q, n)T _1 . Applying the contraction 
rules ( p9|) and then performing the zero-mode integra- 
tion, we reduce the r.h.s. of Eq. (|62|) to the form 

A f 

— J {dp l )dn' 1 w{n ll n' l )(dp 2 )dn' 2 w{n 2l n' 2 ) 

xG R {p l+ )G A (p l _)G R {p 2+ )G A {p 2 _)V{q; n' 2 , n[) 
A 



nr]WQ 



J dn»dn[dn»dn>A(q;ni>,n>) 



xP(q;n , 2 ,n' 1 )A(qr;n' 1 ,n' 1 / ) 



(63) 



Note that in order to simplify this expression, we have 
used time reversal invariance of the classical motion, 



T>(— q; n, n') = T>(q; —n', —n). Substituting ( |63| ) in (54), 
we find the contribution of the (G11G22) term to the wave 
function correlator, 

V 2 C(r 1: r 2 )\ {GllG22) -l 

dri\dri 2 (ni | AT>A\n 2 ) 



{■KV) 2 WI 
1 



dnidn 2 



x{n 1 \AB- 1 A + AB^AB^A 



n 2 



(64) 



The r.h.s. of Eq. ( |64| ) is the sum of the ladder dia- 
grams (corresponding to 1,2,3, . . . intermediate scatter- 
ing processes) yielding precisely the classical propaga- 
tor ns(ri,r 2 ). We see, however, that the first term in 



this sum (corresponding to a motion without intermedi- 
ate scattering) is absent. This term is equal to 

r^— 5- / dnxdn 2 A(ni,n 2 ) 
2 1 (dp)G R (p + )G A (P-) , 



1 



or, in coordinate space 



1 e -\rx-r 2 \/l 

G R (n - r 2 )GA(r 2 - n.) = 



2(ttv) 



■Kp F \rx - r 2 \ 
fo(\ri-r 2 \) 



Including now the leading contribution of the (Gi 2 G 2 i) 
term in (|5^), we get the wave function correlator up to 
the terms linear in the classical propagator /oG^i — r 2 \) 
or in the quantum propagator k q {r\ — r 2 ) = J^(pp\ri 
r 2 \)e 



-\ri—r 2 \/l 



V 2 C( ril r 2 ) ~ 1 + h(r u r 2 ) + k q (n - r 2 ) 



(65) 



where fi = Tig — f . We see that the double count- 
ing problem does not exist anymore. The quantum (k q ) 
and the classical (/1) contributions perfectly complement 
each other, describing the motion before and after the 
first collision respectively. 

Until now, we did not consider the ballistic version 
of the last term in Eq. (|55|). Such a term originates 
from the first order (in T>) correction to (G\ 2 G 2 i). For 
\ r i ~ r 2\ ^ the corresponding contribution to the 
wave function correlator is much smaller than the term 
coming from (GnG 22 ) (second term in Eq. (|65|)), which 
is of order of n^, and thus can be neglected. However, 
this contribution becomes important at \r\ — r 2 \ ~ Af- 
In particular, for r\ = r 2 , the order-D contribution from 
(G12G21) is found to be equal to that from (GnG 22 ), 
yielding 



V 2 C(r,r) ~2 + 2/x (r,r) 



(66) 



After integration over r Eq. ( pq ) determines non- 
universal correction to the average inverse participation 
ratio 



V ( / dr#(r)|- 



1 



[]ng b + 0(l)). (67) 



In the discussion above we did not take into account 
the boundary scattering which determines the classical 
propagator ns(ri,r 2 ) on the scale of the system size 
R but is irrelevant for the matching of the classical and 
quantum contributions on the short scale I <^ R. In prin- 
ciple, one could avoid introducing the additional smooth 
potential and consider the boundary randomness only. 
By analogy with the above consideration, we expect that 
in this case the classical propagator /1 entering Eq. ( |65| ) 
will describe the motion starting from the first collision 
with the boundary. 
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VI. GENERALIZATION TO MIXED BOUNDARY 
CONDITION 

In the case of the mixed boundary condition (Q) the 
trace of the resolvent Tr(K — acquires cuts in ad- 

dition to simple poles which were present for purely dif- 
fusive scattering, a = 1 (see the discussion below). For 
this reason, we cannot write the spectral function S(u>) 
in the Altshuler-Shklovskii form ( j22| ) but have to use a 
more general representation, 



d_ 2 

d 



— ^Trln^ - iuj) 



d(iuj) J 
d 



dte iult Tre- Kt 



d(iuj) 



J dte luJt J drdng{r,n;r,n;t) , (68) 



where g(ri,ni;r2,Ti2;t) is the time- dependent Green's 
function of the operator K characterizing the probabil- 
ity of propagation from the point (1*2,712) of the phase 
space to the point [r\, n\) in a time t. 

Further transformations are straightforward though 
somewhat lengthy. The trace of the Green's function 
in Eq. J68| ) can be written as a sum of the terms with 
n = 2,3, .. . boundary scattering events. Cha nging the 
variables from (r,n) to (9,9',x) (see Sec. [IB ) and per- 
forming the integration over x, we can present the result 
in the form 



Tr(X - iuS) 



2R 

v F 

OO 





9-9' 


/ d9d& 


sin 

2 



(69) 



n=2 

where $ is the integral operator, 

/■2jt 



d9 2 d9' 2 <S>(9 u 9' i; 9 2 ,9' 2 )f(9 2 ,9' 2 ), 







with the kernel 

<Z>(9 1 ,9' 1 ;9 2 ,9 2 ) = 5(9 1 - 9 2 )eiq> 
x {{l-a)6(6x- 9[-9 2 + 9' 2 ) - 



2iujR 
vf 



sm ■ 



(70) 



Physically, $ characterizes the probability of the scatter- 
ing process when a particle moving along the segment 
(6> 2 — > 9' 2 ) is reflected into the segment (£q — > 9[). Rc- 
summing the series ( |69| ) and using that, according to Eq. 



d<t>(9 ll 9' 1 ;9 2 ,9> 2 ] 



2R 

VF 



sm ■ 



$(9u9 



i ! "2, 9' 2 ) , 

(71) 



we rewrite Eq. (|69j) in a compact form 

Tr \n(K - iu) = Tr ln(l - $) , (72) 

up to an irrelevant additive constant (which we drop be- 
low). 

Eigcnfunctions of $ have the form 



(73) 



f(9,9') = e ad g(9-9') , 

so we can write 

Trln(A"- iu) = ^Trm(l- $ z ) , (74) 
l 

where is obtained by restricting the operator <t to the 
space spanned by functions (f73|) with a particular I. The 
kernel of <&/ is given by 



(75) 



9') = exp liW + 2i—\ sin(072)| 

x {(l-a)5(9-9') + ~\sm(9'/2)\} 



where 6 = 9% — 9[, 9' = 9 2 - 9' 2 . We further represent 
as a sum of the contributions corresponding to the two 
terms in the curly brackets in Eq. (|75|), 



$ z (0,0') = <P s l pcc (9')S{9 - 9') + $f iff ((9') . (76) 



Here & spoc is associated with the processes of specular 
reflection, while &f lS corresponds to the diffusive scat- 
tering processes. Using Eq. (|76|), we expand Trln(l — <j>;) 
in powers of $>f lS and then re-sum the series, using the 



structure of the two terms in Eq. ([76|) (the first one is 
proportional to 5(9 — 9'), while the second one indepen- 
dent of 9') to get 



Trln(l-$ ; ) = 

f-271 



2rr 



d91n{l-<f> spec (9)) 



In 



1 



d9(l-<Z> s i pcc {9))- 1 <S>f s (9) 



(77) 



Combining Eqs. @, @, and fn]) we finally obtain 
the following representation for the spectral determinant 
of the Liouville operator, 



Tr hi(K — iu>) = 



1 



d6 



(78) 



x In 



In 



1-- / d9 



1 - (1 — a } e M^R/vF)sin8+2iW 

sin# 



3 — 2i(uR/vp) sin 9 — 2U9 1 -|- , 



The first term in Eq. (|7q) originates from the first 
term in the r.h.s. of Eq. ([77[) and is determined only 
by <I>p cc and not by $f . It thus knows only about the 
motion in a clean system (without boundary scattering) 
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and about the total probability a to be scattered away, 
but not about the differential scattering probability. In 
other words, this term would describe the correlations 
of the density of states if the electrons simply disappear 
(get absorbed) at the boundary with probability a, other- 
wise being reflected specularly. It characterizes thus the 
spectrum of a clean circle (with energy levels broadened 
due to the absorbing boundary). Due to these correla- 
tions, the disorder averaged density of states {v(e)) at 
given energy e fluctuates. These fluctuations are not of 
interest here; one can get rid of them by subtracting the 
(energy-dependent) disconnected part of the level corre- 
lation function, thus modifying the definition of R,2(u>), 

R$(u) = (VA) 2 [{v(e + cj)v(e)) - (u(e + u)){u(e))} , 

(79) 

see Ref. |l8| for an extensive discussion. We mention also 
that, in fact, the first term in Eq. ( |78| ) does not yield 
these correlations fully correctly, since the cr-model is 
not appropriate for .the description of integrable systems. 
Agam and FishmanE3 used Berry- Tabor trace formula to 
calculate the analogous contribution of trajectories not 
scattered by impurities in an integrable system with bulk 
disorder. We do not enter a more detailed discussion here 
since we are not interested in this contribution anyway. 

The second term in Eq. (|7^) describes the disorder- 
induced correlations. Before turning to the calculation 
of the level correlation function, we analyze the struc- 
ture of singularities of the resolvent Tr(K — iuj)^ 1 . As in 
the a = 1 case, it has simple poles lying in the half-plane 
Im lo < which positions £ = itoR/vp are determined by 
the equation 



° / de- 

2 7o < 



sin f 



-2ii0-2£ sin( 



(1 



= 1 



(80) 



However, in contrast to the a — 1 case, the resolvent ad- 
ditionally now has branch cuts induced by zero values of 
the denominator in Eq. (^). These branch cuts can be 
parameterized as 



i 



2sin( 



In- 



1 



2nik - 2iW 



11) 



where l,k are integers, and 9 runs from to n. Physi- 
cally, these additional singularities correspond to motion 
along periodic orbits of the underlying integrable system 
(circle), with the real part Re £ > characterizing the 
total scattering rate out of the orbit. 

To calculate the level correlation function R^oj), we 
use the formulas of Section III, with the spectral func- 



tion S(uj) obtained by substituting the second term of 
Eq. © in Eq. 

Sdoj) = 



2 q2 



R 

vf J dz 



2 hi ■/;(*) | 



Ji(z) = -l + ~ 



dd- 



sm f 



-2il$-2zsin8 _ ft - a) 



(82) 



In particular, the low-frequency behavior has the form 
( pl| ) with the coefficient A given at a <C 1 by 



A = 



256 
9^ 



3.12 



(83) 



As previously, this result is valid until the last term (l/gb 
correction) in Eq. (21) becomes of order unity. This hap- 
pens at a frequency which plays the role of the Thouless 
energy. For a <C 1 this energy scale is E c ~ avp/R, 
which is the inverse of the characteristic relaxation time 
at a <C 1. The a- model only applies when the Thou- 
less energy is much larger than the mean level spacing, 
i.e. for a. S> g^ 1 . For even lower values of a the system 
exhibits integrable behavior. 

As usual, at frequencies above the effective Thouless 
energy E c the level statistics are totally different from 
the RMT predictions. To demonstrate this, we calculate 
the smooth part of the level correlation function for in- 
termediate frequencies E c -C u> <C vf / R using Altshulcr- 
Shklovskii formula (p8|). We first notice that the contri- 
butions of all Si with I ^ are suppressed by the small 
factor uR/vp as compared to the term So, and thus can 
be neglected. Evaluating the integral in Eq. (^2|) at / = 
and extracting the leading-order term, we find 



3a 2 
ATt 2 gl \luR 



(— V 



In 



/ ujR 

\OtVF 



O(l) 



(84) 



The result (|8j) matches Eq. (|8J) at w r 
and drops sharply for higher frequencies. 



E c — avp/R 



VII. SUMMARY AND DISCUSSION 

In this paper we have used the ballistic cr-model ap- 
proach to investigate statistical properties of energy lev- 
els and eigenfunctions of a circular billiard with diffusive 
surface scattering. For this simple model of a chaotic 
system we calculated explicitly non-universal deviations 
of the statistical properties from the random matrix the- 
ory. These non-universal properties are determined by 
the classical dynamics and turn out to be very different 
from the two examples available in the literature - diffu- 
sive systems and ballistic systems with bulk ^-correlated 
disorder. We believe that our results are not specific 
for a particular model considered but rather reflect non- 
universal features of generic chaotic systems and their 
qualitative difference from the corresponding properties 
of systems with bulk shot-range impurities. Below we 
summarize our main findings. 

The spectral statistics (Section |l|) deviate from its 
RMT form on a frequency scale set by the inverse flight 
time (playing the role of the Thouless energy). At higher 
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energies, the level number variance saturates and oscil- 
lates in agreement with predictionso for a generic chaotic 
system. Surprisingly, the two-level correlation function 
shows non-decaying (though weak) oscillations with the 
period of the mean level spacing at high frequencies, pro- 
ducing a t)-like spike in the spectral form-factor at the 
Heisenberg time. In Section VI the analysis of the spec- 
tral statistics is generalized to the case of a mixed bound- 
ary condition when the "Thouless energy" is parametri- 
cally smaller than the inverse time of flight. 

In Section [V we presented a thorough study of the 



parametric level statistics, with magnetic field playing 
the role of the external parameter. We identified all rel- 
evant regions in the frequency-magnetic field plane and 
calculated the parametric two-level correlation function 
and the parametric level number variance in all of them. 
In particular, a surprising result was obtained in the re- 
gion of high magnetic fields, where the parametric corre- 
lation function was found to be independent of the mag- 
netic field. In other words, the density of states retain 
finite memory even in the limit of very large fields after 
levels have undergone arbitrary many avoided crossings. 

In Section |y| we analyzed spatial correlations of eigen- 
function intensities. Since naive application of the bal- 
listic u-model led us to the problem of double counting 
(with one and the same contribution appearing twice - 
classically and quantum- mechanically) , we had to rean- 
alyze the tr-model derivation. For this purpose we intro- 
duced a smooth random potential and obtained the tr- 
model from averaging over it, following Refs. |5^ (see Ap- 
pendix). We have found that while in the limit r/r tr — > 
(with t -1 and Tj" 1 being total and transport scattering 
rates, respectively) the obtained action takes the conven- 
tional form ( [A6| ) of the ballistic tr-model, the behavior of 
the propagator is different at short distances. This af- 
fects the wave function correlations at short spatial scales 
r < vft. The final result has a form (|6^) with the quan- 
tum (k q ) and classical (/i) terms corresponding to the 
motion before and after the first collision respectively. 

In the present context of a system with diffuse bound- 
ary scattering, introduction of an additional smooth ran- 
dom potential satisfying d <C I <C R <C kr can be consid- 
ered as a technical trick allowing to obtain the tr-model 
in the conventional form (Q). In principle, it should be 
possible to derive the tr-model directly by averaging over 
the boundary disorder, but in this case the action will be 
more complicated, since the system size will essentially 
play the role of opr. Let us stress that the additional ran- 
dom potential does not affect the results for the energy 
level statistics and for the smoothed eigenfunction cor- 
relations. On the other hand, the corresponding mean 
free path I = vft manifests itself explicitly in Eq. j65| ) 
for the eigenfunction correlations by setting the scale at 
which the Friedel-type oscillations get smeared. 

We believe that averaging over an additional smooth 
quantum random potential is of conceptual importance 
for the problem of non-universal features in the level 
and eigenfunction statistics of a conventional chaotic bil- 



liard (without boundary scattering). Indeed, a consen- 
sus seems to have been reached by now that the energy 
averaging by itself is insufficient, and one has to aver- 
age over some class of systems (with the same classical 
dynamics) in order to derive the tr-modclEZI. Further- 
more, without such an averaging one cannot detect n<m- 
universal features, since the statistics are not sufficientEa. 
A smooth quantum random potential with parameters 
chosen in such a way that VFTtr/R oc hT a with a > 
and vft/R oc h with < b < 1 is exactly the required 
type of ensemble averaging. With this averaging, the 
ballistic tr-model can be rigorously derived. Let us em- 
phasize that the ballistic action obtained in this way will 
have the conventional (obtained by gradient expansion) 
form (Q) only at length scales exceeding I = vpr. If one is 
interested in eigenfunction correlations at shorter scales, 
one should avo id t he gra dient expansion and use the more 
general form (A5), (A8), as explained in Sec. [v] and in 
Appendix. 

An additional smooth random potential discussed here 
resembles the one introduced by Aleiner and Larkin in 
Refs. ^,^6| where the problem of weak localization in bal- 
listic chaotic systems was considered. There is, however, 
a conceptual difference: while Aleiner and Larkin intro- 
duced fictitious disorder in order to mimic diffraction on 
boundaries of the billiard, we consider a real random po- 
tential, i.e. an ensemble of systems, averaging over which 
allows one to derive the tr-model, as explained above. 

We close the article by mentioning two open issues. 

(i) The problem of repetitions! (which appear not to 
be counted properly in the tr-model approach) still awaits 
its resolution. For the model considered in the present 
article this problem does not apply since in view of the 
diffusive nature of boundary scattering, all directions of 
motion after a scattering event are allowed so that the 
repetitions are irrelevant. 

(ii) We used the boundary condition for diffusive scat- 
tering in a linear form, i.e. we supplemented the Liouville 
operator determining the quadratic form of the action 
by the boundary condition. This was sufficient for the 
problems considered in the present article, when the rele- 
vant tr-model correlation functions are determined by the 
structure of the action in the vicinity of spatially homo- 
geneous configurations and therefore the results are gov- 
erned by the eigenfunctions and eigenvalues of the Liou- 
ville operator. However, in general, this is not sufficient, 
and a boundary condition on the g(r, n) field is needed. 
In particular, one would need such a general boundary 
condition to calculate "tails" of various distribution func- 
tions (of relaxation times, eigenfunction amplitudes, local 
density of states, inverse participation ratio, etc.), analo- 
gously to how it has been done for diffusive systemsE3Q. 
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APPENDIX A: BALLISTIC cr-MODEL FROM 
SMALL- ANGLE SCATTERING 

Here we outline derivation of the cx-model based on av- 
eraging over a smooth random potential U(r) (see Sec- 
tion |yp, following Ref. |55|. After the averaging and the 
Hubbard-Stratonovich decoupling by a supermatrix field 
Q(r,r') one finds the action 

F[Q] = Str hx[E + i^A - H - Q(r,r')V(r - r')] 



where w(n,n') = 2ttvV(pf \n — n'\) is the scattering 
cross-section. Performing now the gradient expansion of 
Eq. (A5) and using r tr 3> r, one gets the action of the 
ballistic cr-model (generalizing Eq. (|l|) to the case of a 
non-isotropic disorder scattering), 

F = irvv F j drdn Str AT _1 (r, n)nVT(r, n) 

+ ~i~ J drdndn ' Str Q{r,n)w(n, n')Q(r,n r ) . (A6) 



Let us emphasize that the action (A5) takes the form 
(|A6| ) only in the limit r/r tr — > 0. While for most pur- 
poses the difference between these two formulas is irrel- 
evant, it is of crucial importance for settling the double 
counting problem considered in Sec. since it is related 
to the short-scale behavior of the cr-model propagator. 
For this reason we use there the action in the form (A.5) 
(and not the approximation ( |A6| )). 

To illustrate the con nection and the difference be- 
tween the actions (A5) and (|A6|), it is instructive to 
write down their quadratic forms at low momenta. 
We write T(r,n) = 1 — W(r,n)/2 and introduce 
the angular harmonics of the VF-field, W(r,m) = 
J dnexp(— im<fi n )W (r , n), and of the scattering cross- 
section, w m — f (d<fi/2ir)e" n ^w(d). The resailt for the 



X r section, w m — J [dcp/ Zn)e ^wup). ihe resaiJ 

+ -Str / drdr'Q(r, r')V{r - r')Q(r', r) , (Al) quadratic form of the action © then readsEl 



where Hq = p 2 /2m is the free Hamiltonian. The corre- 
sponding saddle-point equation has the form of the self- 
consistent Born approximation (SCBA) and possesses 
a set of translationally-invariant solutions <2(r, r') — 
Q(r — r'), which can be most conveniently written down 
in the momentum space, 

Q(p) = ReG R ( P ) + iTAT-HmGnip) . (A2) 

Here Gr(p) is the retarded Green's function in SCBA, 



= \ j {dq)StxW^ l2 {-q)T rnml {q)W mlai {q) , (A7) 



where 

r mm '(q) - 

iv F w m w m 



— (too - w m )5 mm > + - — q S m0 5 m > 
w 2w 



(qS-n 



m,m' ' + 1 j 



(A8) 



2tv 



G R (p) = IE - H (p) 

i = - / (dpJVip-pJlmGuipx) , 
2t p J 



(A3) 



where q — q x + iq v . On the other hand, the quadratic 
terms in the action (A6) are given by Eq. (A7) with the 
kernel T replaced by 



and the matrices T belong to the super coset space 
[/(2|l,l)/f/(l|l) x C7(l|l). Allowing for slow variation 
of T with the spatial coordinate r + = (r + r')/2 and 
with the direction n of the momentum p yields the soft 
modes, 



fmm'(g) = 7T*4(w - Wm)^ 



(A9) 



Inverting Eqs. (A8), (A9) at small c/, one finds that the 



Q(r+,n) =T(r + ,n)AT- 1 (r+,n) 



(A4) 



corresponding propagators T> = T 1 and D = T 1 are 
identical, up to a constant term, 



with Q 2 (r + ,n) = 1. The action for these soft modes has 
the form (we set u> = 0) 



F = Sir Id { E - - H a - - / dnQ(r,n')w(n,n'] 



drdndn' Str Q(r,n)w(n,n')Q(r,n') , (A5) 



V mm >(q) - V mm >{q) 



1 



(A10) 



The physical meaning of this difference becomes clear 
from the calculation in Section [v| (avoiding the momen- 
tum expansion). Specifically, the propagator for the ac- 
tion S is given by a series of ladder diagrams, beginning 
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from the term "with —1 scattering" [the term B 1 in 
Eq. (f30|)], while that for the action S starts from the term 
with zero collisions (free motion), i.e. from the second 
term on the r.h.s. of Eq. @. As a result, V = V-B~ 1 , 
which reproduces exactly Eq. ( |A10| ). 
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